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CANCEROUS TUMOR DETECTION USING OPTICAL 
VASCULAR FUNCTION 

This invention was made with Government support under Grant DAMD 17-02- 1-0570 
awarded by the U.S. ARMY MEDICAL RSCH. ACQUISITION ACTIVITY . The 
Government has certain rights in this invention. 



INTRODUCTION AND HYPOTHESES 

The long-term goals of this research are to develop inexpensive technology to 
Improve sensitivity and specificity (lower false-negative and false-positive rates) for early 
breast cancer detection and diagnosis. We believe that enhanced functional 
(physiological) optical imaging using a new type of contrast based on the unusual 
vascular function of tumors (atypical oxygenation improvement, atypical vaspactivity. 
and blood pooling) can meet these goals. Improved imaging through dense breasts 
(where X-ray mammography is less successful) is a secondary goal of this work. To our 
knowledge we are the only group investigating differential optical imaging using the 
vasoactivity of oxygen (O2) and carbon dioxide (CO2). We have been investigating this 
differential vasoactive optical imaging (DVOI) approach in animal model studies 
supported by an IDEA grant from the U.S. Army Breast Cancer Research Program. 
That work has demonstrated strong contrast between cancerous and noncancerous 
tissue during differential imaging in rodents in association with inhalation of O2/CO2 gas 
mixtures, providing proof-of-principle preliminary data in support of this TRC award 
application. 

The contrast achieved by DVOI results from the vasculature in tumors and can arise 
from atypical oxygenation improvement, atypical vasoactivity, and blood pooling, as 
monitored by varying the levels of inspired O2 and CO2. We will use these differential 
vascular function measurements to augment the cancer-specific static contrast derived 
from elevated hemoglobin concentrations from angiogenesis and reduced local 
hemoglobin oxygenation from tumor hypoxia. A single imaging DVOI system can 
monitor both static and dynamic contrast mechanisms, thus providing the best possible 
sensitivity and specificity from an optical imaging system. CO2 and O2 are attractive 
contrast-enhancing agents because they are safe and require no injection or lengthy 
times between administration and imaging. 

We hypothesize that the specificity of the differential contrast available with the DVOI 
approach will be sufficiently great to allow tumor detection with higher sensitivity, even at 
the poor spatial resolution available using optical imaging through the human breast. 
Other imaging techniques such as functional MRI (fMRI) or positron emission 
tomography (PET) are powerful imaging modalities, even with their low spatial 
resolutions. The strong contrast observed in our animal model studies and the finding 
that optical techniques can detect and locate picomole variations in chromophore 
concentrations over optical thicknesses comparable to those of the human breast^ 
further supports this hypothesis. 

SPECIFIC AIMS 

We have proven the capabilities pf DVOI for detecting cancer in an animal model. 
We propose the translation of this work into the clinic to determine whether the DVOI 
approach is equally useful in humans. Our specific project aims are to: 
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(1) Modify and improve our DVOI system to allow human imaging. 

(2) Perform clinical studies of DVOI mammography imaging on women scheduled 
for biopsies as a result of screening mammograms or suspicious lumps. 

(3) Perform limited imaging on animal models to understand variability in the 
vasoactive response better and to test refinements in the imaging system, 
protocols, and data analysis. 

(4) Modify and refine gas inhalation protocols for women by optimizing variations in 
gas composition and timing of gas inhalation to determine whether differential 
contrast can be further enhanced. 

(5) Implement data analysis procedures to facilitate identification of significant image 
regions, to assist in understanding the large amount of data obtained in time- 
dependent two-dimensional (2-D) imaging, and to identify possible lesions for 
breast imaging. 

BACKGROUND AND SIGNIFICANCE 

Significance 

Our objective is to test the capability of DVOI to meet the strong need for improved 
noninvasive imaging techniques for breast cancer detection. Improved imaging through 
dense breasts is a secondary goal of this work. Detecting cancer eariier and reducing 
false-negative results will improve clinical outcomes, and lowering the number of false- 
positive results will decrease the cost and physical and emotional trauma associated 
with unnecessary biopsies. X-ray mammography, the standard screening modality for 
breast cancer, is less effective at detecting cancer in younger women's breasts, which 
are denser than those of older women. Although the risk of carcinogenesis resulting 
from X-ray mammography is relatively low, concerns about risks of exposure over many 
years of screening are valid. For these reasons, other imaging techniques are being 
used and studied to augment X-ray mammography, including ultrasound, MRI, Tc-99m 
sestamibi scintimammography, optical techniques, and PET. 

Optical techniques offer advantages for breast imaging, including functional imaging 
(imaging that provides information on tissue state and function), inexpensive 
instrumentation (our current imaging system costs less than $1,000 for parts), and no 
use of ionizing radiation. Optical imaging could prove useful as a secondary imaging 
modality to X-ray imaging for diagnosing, staging, or monitoring treatment of breast 
cancer If DVOI proves successful, it could also find use as a primary screening 
modality. Because of its simplicity and low cost, DVOI can be efficiently incorporated 
into an X-ray or ultrasound imaging system^ to provide functional information to 
complement the physical imaging of these modalities. DVOI may prove more effective in 
imaging dense breasts and mav reduce or avoid the sometimes painful compression 
used for X-ray mammography.^ 

Optical Breast Imaging 

Optical imaging through the human body was one of the first modalities used for 
medical imaging, with early studies performed in the 19*^ century. ^'^ Optical imaging of 
the breast was reported in 1929.''° Optical mammography was closely studied in the 
1970s and 1980s.^^'^^ although it proved inferior to X-ray mammography.^"*-^^ The 
primary problem with optical mammography is spatial resolution. Optical mammography 
has a spatial resolution of 0.5 to 1 cm, which means that blurring reduces contrast in 
smaller tumors. This limitation can be overcome by providing functional imaging 
information. 
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Interest in optical techniques for detecting breast cancer has increased » spurred by 
better understanding of liaht propagation through tissue,^®-^^ time-resolved optical 
techniques for imaging^'^' and quantifying tissue properties,^^'^^ and new image 
reconstruction algorithms.^^^ Industry (e.g., Carl Zeiss,^® Phillips,^ Seimens.^) has 
developed prototype imaging systems Applications of advanced optical mammography 
techniques have shown promising results.^ '*'^'^®'^^ A number of groups are examining 
the use of contrast agents for breast cancer detection. ^^"^ Other useful optical imaging 
studies have incorporated window chambers. ^^'^^ A number of ex vivo studies of the 
optical properties of normal and diseased breast tissue have been performed.^® Our 
studies complement these studies by investigating a new and promising method of 
increasing contrast. 

Functional Optical Imaging 

Conventional imaging/static contrast. Whereas X-ray imaging primarily provides 
structural information, optical spectroscopy imaging can provide infonmation both on 
structure and tissue function. For example, optical measurements at different 
wavelengths can indicate total hemoglobin content and oxygenation — afunctional 
information that is significant for breast cancer detection. Tumor angiogenesis typically 
leads to elevated local hemoglobin concentrations.^®'^® In addition, tumors are often 
hypoxic, which can be observed optically as a decrease in hemoglobin oxygenation. 
Because tumors that are more hypoxic tend to be resistant to radiotherapy and 
chemotherapy and are more likely to be metastatic or invasive.^^"^^ the degree of tumor 
hypoxia can be used to guide treatment. Tumor morphology also provides a source of 
contrast through variations in the optical scattering coefficient.^ We propose to use the 
full range of available optical contrast by augmenting functional optical imaging using 
differential measurements related to tumor vascular function. We believe the broadest 
use of available contrast will be most effective for improving sensitivity and specificity. 

Differential vasoactive contrast The atypical characteristics of vasculature 
produced through tumor angiogenesis^^ provide the basis for our DVOI approach. 
Blood vessels in tumors often exhibit distended capillaries with leaky walls and sluggish 
flow. These properties provide at least three types of contrast for optical imaging in 
conjunction with varying levels of inspired O2 and CO2. These types of contrasts are due 
to atypical oxygenation improvement, atypical vasoactivity. and blood pooling ."^^^ Dr. 
Britton Chance and collaborators have observed further evidence of compromised 
microvascular control through measurements of NADH (nicotinamide adenine 
dinucleotide) signals.^ Because both O2 and CO2 are vasoactive, atypical tumor 
vasoactivity arising from administration of changing levels of these gases should provide 
strong imaging contrast. Tumor vessels are often contorted and leaky; thus, blood 
pooling in these vessels will delay response to oxygenation changes, providing another 
good contrast mechanism. Blood pooling itself can contribute to the atypical 
oxygenation improvement in tumors. However, as our preliminary studies have 
indicated, atypical oxygenation improvement persists beyond the transient response 
caused by blood pooling. 

The unusual vasculature in tumors can be readily measured using functional optical 
imaging. For example, opposing vasodilation and vasoconstriction responses after 
15% CO2 and 85% O2 (carbogen) inspiration are readily detectable by comparing 
hemoglobin content before and after carbogen is administered. Similarly, the changing 
response in tumor oxygenation after increased O2 administration is easily measured by 
monitoring hemoglobin oxygenation levels before and after the O2 level is increased. 
Changes associated with blood pooling are observable in delayed oxygenation changes 
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in the tumor. DVOI could also be combined with quantitative measurements of oxy- and 
deoxy-hemoglobin to Improve overall sensitivity and specificity. 

DVOI has a good potential for providing functional discrimination between benign 
and malignant lesions. Benign lesions tend to have rounded vasculature while 
malignant lesions tend to be more angular. Because the vasculature is different, it is 
likely that the vascular response to Oa and COz will also be different. If DVOI can 
provide information distinguishing benign and malignant lesions, this would be valuable. 

There are additional motivations for examining differential contrast such as that 
associated with tumor vascular function. First, because the breast is highly 
heterogeneous, comprising the lobes (glandular tissue), fat, connective tissue, ducts, 
and supporting vasculature, using a broader palette of contrast mechanisms should 
provide more specificity for optical imaging and help compensate for that heterogeneity. 
Second, the more successful noninvasive optical measurements (e.g., pulse oximetry, 
functional brain imaging) are differential or dynamic. In pulse oximetry the 
differential signal of the cardiac pulse is combined with the optical signals of oxygenated 
and deoxygenated hemoglobin to derive blood oxygenation levels without detailed 
knowledge of the tissue type and level of scattering. For brain imaging,^^ differential 
measurements of hemoglobin or hemoglobin oxygenation are performed while the 
patient performs specific tasks such as finger tapping or is visual stimulated. As a result 
of the successes of these approaches, we will also combine dynamic and functional 
measurements for the studies proposed here. Finally, recent theoretical work has 
demonstrated improved results using dynamic or differential optical imaging techniques, 
both of which rely on changes in optical contrast over time. 

PRELIMINARY STUDIES 

Differential Vasoactive Imaging Studies 

With support from the U.S. Army Breast Cancer Research Program, we are studying 
DVOI for mice and rats . Animal models allow us to monitor contrast for a range of 
tumor sizes and stages of development. Furthermore, the measurements are 
noninvasive and thus can be readily repeated on animals as our instrumentation and 
methods are refined. We perform measurements through tissue thicknesses similar to 
those of the human breast by partially immersing the anesthetized animals in liquid 
tissue phantoms that simulate the optical properties of breast tissue. Although this 
approach does not allow for the effects of tissue heterogeneity in the breast, it is the 
most practical method for studying contrast before beginning studies with human 
subjects. 

We have constructed a continuous wave (CW) immersion imaging system to perform 
DVOI. Figure 1 provides a schematic of the instrument used for immersion imaging. 
The imaging system is composed of a near-infrared (NIR) light source, an immersion 
box, and a camera. The light source is made up of an array of bright light emitting 
diodes (LEDs) that emit NIR radiation with peak intensities at either 780 or 840 nm 
(Epitex L780-01AU and Epitex 840-01 KSB, respectively). Switching between LED 
arrays enables measurements at different wavelengths and the determination of 
hemoglobin content and hemoglobin oxygenation. This light source is directed at the 
sample immersion box, which contains the study anihial in a heated (S/^'C), matching 
medium composed of water, ink, and Ropaque (Rohm and Haas Company). (Ropaque 
is made up of submicrometer polymer spheres used as an alternative to titanium dioxide 
in paints.) This immersion medium approximates the scattering and absorptive 
properties of the mouse tissue. The front of the immersion box is imaged onto the 
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camera. Images at each individual wavelength are then collected, digitized (8-bit 
resolution), and sent to the computer for analysis. 

The compensation provided by immersing the animal in a tissue phantom improves 
image quality by removing changes in contrast associated with changes in tissue 
thickness and geometry, allowing better use of the dynamic range of the camera and 
providing more uniform illumination. When the match is good, the tissue almost 
disappears, and the image shows variations due to internal structure and contrast, which 
is what we want for in vivo imaging. The immersion medium serves to: (1 ) allow study 
of an effective tissue as thick as is typical for the human breast, and (2) enhance our 
measurements by eliminating the effects of boundaries. Although the tissue phantom 
lacks the heterogeneity of the human breast, there is considerable heterogeneity in the 
animal itself. Nonetheless, the effects of breast heterogeneity and human vasculature 
can be examined only in human studies, as proposed below. 

Inrmension 
Box 



0 




Corrputer 



Figure 1. Schematic diagram of the CW Immersion imaging system. 

Tissue phantoms are prepared using our established methods. ^'^^ An initial tissue 
phantom is prepared, an animal is immersed between the source and collection fibers, 
the changes in amplitude and phase are measured, and the phantom composition is 
adjusted according to the optical properties determined from the immersion 
measurement. This process is repeated until the optical properties of the immersion 
medium and the imaged tissue agree to within a few percent. The thickness of the 
tissue phantoms is varied by inserting Plexiglas sheets into the box containing the tissue 
phantom for the CW measurements. 

We have compared three cameras for the CW imaging: a Dragonfly CCD (charge- 
coupled device) camera (Point Grey Research), a Pulnix TM-9701 CCD camera coupled 
to a Stanford Photonics Gen III image intensifier, and an ImagingSource DMK-3002-IR. 
The digital Dragonfly camera offers a significant improvement in signal-to-noise ratio 
(S/N) over other video cameras. Although the Dragonfly has a lower absolute sensitivity 
in the NIR region compared with the other video cameras, it has lower read noise and is 
capable of longer exposure times (>60 s), which will be important for imaging thicker 
tissue samples. We increased light throughput onto the imaging sensor by 20% by 
installing a large-aperture lens with high NIR transmission (JML Optics). 

Animal Models 

Human breast cancer cells (MDA 231) and mouse embryonic fibrosarcomas were 
grown in Dulbecco's minimum essential medium (DMEM) with glutamine and 10% fetal 
bovine serum. The cells were harvested when they were 80% confluent, using 0.25% 
trypsin. Cells were injected subcutaneously on the dorsum of the female athymic nude 
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mice (approximately 23 g. Harlan Laboratories). Both cell lines were used at a 
concentration of 2-3 million cells in 100 \i\ of DMEM for each animal. The tumor 
volumes were measured twice weekly. 

Animal Imaging 

Imaging experiments were conducted on animals with tumor volumes of 500-1000 
mml We used two-four animals for each experiment. After being anesthetized with 40 
mg/kg of pentobarbital, the mice were secured to a 3-mm Plexiglas platform with black 
vinyl tape. Anesthesia was given in further doses of 20 mg/kg as needed to reduce 
stress associated with immersion and to keep the animal immobilized. Carbogen or air 
was administered to the immersed mouse via a nose cone at a flow rate of 
approximately 3 l/min. The optical path length of the immersion box was adjusted to 
match the thickness of the mouse ('<-2-2.5 cm). At this thickness, the exposure time of 
the camera allowed us to measure both wavelengths at approximately three frames per 
second. At the end of the experiment, the animals were sacrificed with an overdose of 
anesthetic agent. 

Images of individual mice were recorded before, during, and after the administration 
of carbogen. Figure 2a shows one of these static images taken 134 s following the 
administration of the carbogen. The approximate outlines of both the mouse and the 
tumor have been placed on top of the image as a guide. The mouse's head is out of the 
immersion medium and is above the field of view in Figure 2. The hind legs and tail are 
seen at the bottom of the image. Figure 2b shows this same image after the subtraction 
of a background, which is simply an image of the mouse before the carbogen was turned 
on. Although the boundaries of the mouse and tumor are obscured by the good match 
with the immersion medium, it is clear from the difference image in Figure 2b that there 
are distinct regions of contrast between the tumor and the surrounding tissues of the 
mouse. 




Figure 2. Images of mouse at 840 nm at 134 s after administration of 
carbogen; 

(a) static image; (b) image from 2a with background subtracted. 

Temporal Variation in Differential Contrast 

The enhanced contrast between the tumor tissue and the mouse tissue due to the 
inhalation of the carbogen was monitored by averaging the changes in intensity over 
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areas within the difference images. Figures 3 and 4 show these averaged data for 
differences in the 780 nm and 840 nm images, respectively. The squares represent 
changes in the tumor tissue, the circles indicate an adjacent region within the mouse that 
does not contain the tumor, and the line represents the average of a part of the image 
not containing the mouse. The maximum change for both wavelengths is approximately 
±10 units, and it is clear from the figures that distinct differences occur for the dynamics 
of the tumor tissue when compared with the normal mouse tissue. Furthermore, the 
background, which is a measure of lower limits for detection, varies just ±0.2 units. 

Figures 3 and 4 indicate that several regions (e.g., near 55 s at 780 nm, and near 
135 s at 840 nm) show strong contrast between tumor and surrounding tissue. 
Additional contrast is found after the carbogen is stopped; for 840 nm, the relative 
intensity of tumor and surrounding tissue reverses. 
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Figure 3. Temporal evolution of regions of the difference 
images at 780 nm. 



Figure 4. Temporal evolution of regions of the differ 
images at 840 nm 



Although images at a single wavelength such as Figure 2b can be useful for cancer 
detection, it is also of interest to determine the changes in oxyhemoglobin and 
deoxy hemoglobin. We have analyzed the same image data set used to produce Figures 
3 and 4 to calculate approximate path-integrated oxyhemoglobin and deoxyhemoglobin. 
The absorption at 780 nm and 840 nm can be described as: 

where A is the wavelength of interest, [Hb] and [HbOJ are the concentrations (moles/L) 
of deoxygenated and oxygenated hemoglobin, respectively, and ^ is the molar 
absorption coefficient.^ Using Beer's Law, we can describe the change in the 
absorption coefficient at time t after a baseline image has been taken as: 



A/i,"=//f-//,'**«^"'' = 2.31og, 



' baseline 



II 



where / is the intensity of transmitted light and / is the pathlength in cm, corrected 
appropriately for the differential pathlength factor for the animal tissue.^ ®^ We can 
obtain a rough measure of the change in path-integrated oxyhemoglobin and 
deoxyhemoglobin concentrations by assuming that the differential pathlength factor is 
the same at both wavelengths. By manipulating equations (1) and (2), we see that: 
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23 
/ 



Because of the finite bandwidth of the LEDs. we calculated the absorption coefficient 
by integrating the wavelength-dependent absorption coefficient with the normalized 
spectra of the LEDs for each wavelength respectively; 



This led to the following equations for the concentrations of Hb, Hb02, Hbtotat at time t: 



A[Hbtotlt) = A[Hb]it) 4- A[Hb0^lt) 

We used these calculations to calculate the approximate temporal variation of the 
total hemoglobin, oxyhemoglobin, and deoxyhemoglobin shown in Figure 5. These 
values were in turn used to calculate the approximate change in O2 content 
(oxyhemoglobin change, minus deoxyhemoglobin change) shown in Error! Reference 
source not found.. Several observations arise from these images: The tumor 
vasculature shows more erratic behavior, as seen from the oscillations at the beginning 
of carbogen inhalation, the failure to return to baseline for the total hemoglobin 
concentration (Figure 5), and the overshoot in O2 content at the end of the carbogen 
inhalation (Errorl Reference source not found.). The magnitude in changes of 
oxyhemoglobin and deoxyhemoglobin are accentuated in the tumor (Figure 5). The 
increase in O2 content of the tumor is delayed relative to the rest of the animal (Figure 5 
middle and Errorl Reference source not found.), which may be due to blood pooling in 
the tumor. 





4HbO,lt)^23^ - 5.255 MO-^*log,o-^ + 7.996* 10-^ Mog.0-^ //, 
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Figure 5. Temporal variation of relative changes in 
total hemoglobin (top), oxyhemoglobin (middle), and 
deoxyhemoglobin (bottom) during carbogen inhalation. 

The tumor region is shown by the dashed line; the 
region on the mouse torso away from tumor is shown 
by the solid line. 




Figure 6. Temporal variation of relative changes in total O2 content (oxyhemoglobin change, 
minus deoxyhemoglobin change) during carbogen inhalation. The tumor region is shown by the 
dashed line; the region on the mouse torso away from tumor is shown by the solid line. 
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The same processing used for Figures 5 and 6 can be used to produce images 
representing approximate path-integrated oxyhemoglobin and deoxyhemoglobin (Figure 
7). These differential vasoactive images show a dramatic increase in tumor contrast as 
compared with a raw or static image (Figure 2). 





Figure 7. Relative concentrations of oxyhemoglobin (left) and deoxyhemoglobin (right) 
concentrations at 140 s (100 s after carbogen administration). 

Principal Component Analysis 

The imaging experiments described above generated large sets of data. Typically, 
images with 10 pixels at two wavelengths are recorded every 2-10 seconds over the 
cycling period of carbogen administration (approximately 10 to 20 minutes). From our 
preliminary studies with this method, we expect to see <7% change in image intensity 
following carbogen administration. Because extracting such small signal changes from 
large data sets poses a formidable challenge, researchers have developed techniques 
that generate smaller sets of orthogonal images to describe the generated data.®^'^^ In 
practice, these methods have been shown to accurately describe data sets of 10,000 
images with only --1 00 eigen images.®^ 

In the most basic adaptation of these methods, known as principal component 
analysis (PCA), the set of recorded images is represented by. ^'^^ 

/ = /(^x) 
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where x describes the spatial pixel grayscale values of the image, and t is the time at 
which the image data was collected. Researchers have shown that these images, f(tx), 
can be decomposed into the set of orthogonal functions an(t) and ^fxj by: 

A series of Ttime images containing P pixels can be described by the matrix: 

V(U) /(U) ... f{\.p)' 

/(2.1) /(2,2) ... /(2,P) 



M = 



... f{T.P) 



This matrix can then be decomposed into the different art(t) and (pnlx) components 
through the general technique of singular value decomposition: 



,and U = 



and 



M=AUV^ 



The columns of V contain the orthonomnal spatial basis functions, the orthonomial 
columns of A describe the time-dependence of the spatial basis functions, and U 
contains the weighting factors for the two matrixes A and V. We have applied this 
simplified PCA method to our data as a first step toward simplifying the extraction of 
spatial images. ' 



As a first step in processing the data, we 
followed this procedure to determine changes 
in oxyhemoglobin and deoxyhemoglobin. 
scaled by some pathlength factor / as described 
above. The time-dependent images that 
describe d[Hb] and ^[HbO^ were ordered into 
a matrix as shown in equation (10). and the 
singular value decomposition was carried out to 
obtain the matrices A, U, and V. Figure 8 
presents a plot of the normalized scaling factors 
contained along the diagonal of U. Only the 
first three or four eigen images contribute 
significantly to the set of images that describe 
the hemoglobin dynamics in our study. 



NonnaKged Eigeiwahies 




Figure 8. Normalized eigen value spectrum. 



Figure 9 shows the first two eigen images corresponding to the first two columns of 
matrix V. The contrast between the tumor and the surrounding tissue is evident in the 
second image. The time-dependent weighting of the second eigen image in the ls[Hb](t) 
and L[HbOii(t) sets of images can be determined from the matrix product of A*U, and is 
shown in Figure 10. 
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Figure 9. First two eigen images from principal component analysis. 
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Figure 10. Temporal variation of the eigen image scaling factor. 
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Remaining Project Goals 

Our U.S. Army Breast Cancer Research Program IDEA award ends in early May 
2004. The primary unfinished goal of this project is the comparison of differential 
contrast with static contrast. In the remaining 5 months of this project, we will be 
performing frequency domain measurements to examine static contrast. These 
measurements will be performed using an existing homodyne frequency domain system. 
Source and detector fibers will be scanned in the immersion medium to enable 
quantification of the static contrast and provide semiquantitative results from the 
differential measurements. 

IVIethods for Improving Quantification (possible future work following proposed 
project) 

Dr. Fans has extensive experience with the frequency domain technique,^*^'®^^^ a 
method used to improve the quantification of tissue properties, and has more than 1 5 
years of experience with optical tomography instrumentation^^*^"^^^ and algorithms.^^ 
This experience could prove useful in enhancements of this technique at some future 
date, if necessary. 



p. 13 of 28 



Applicant Reference No. US-51 13-2P 



RESEARCH DESIGN AND METHODS 
Aim 1. Enhance Imaging System 

We will modify our existing imaging system to improve its performance by: adding 
another wavelength to enhance the imaging of water, increasing the illumination power, 
and increasing camera sensitivity. These modifications will enable us to image through 
large tissue phantoms with S/Ns limited only by shot noise, which is a fundamental 
limitation for any imaging process. High S/Ns can be very effective for differential 
imaging because image heterogeneity is removed during the image subtraction process. 
That is. subtraction of two images taken of the same field of view yields an image of zero 
intensity if nothing has changed. 

Water concentrations are known to influence measurements of hemoglobin.®' Thus, 
performing imaging at a wavelength dominated by water absorption may assist in 
quantifying oxyhemoglobin and deoxyhemoglobin measurements. Because of the high 
fraction of water in blood, images with dominant water absorption may also help monitor 
blood volume directly. Although the change in water content associated with 
vasodilation or vasoconstriction is relatively small, we have found that the differential 
imaging is quite sensitive to such changes. Thus, we may be able to monitor changes in 
blood volume directly using differential images at 970 nm, a wavelength dominated by 
water absorption. Observations that vasoconstriction may lead to changes in 
hematocrit®^ may mean that a water-based measurement of blood volume can also 
provide information on blood plasma changes, which are somewhat different from the 
changes provided by hemoglobin measurements. Measuring blood plasma changes 
may prove to be infeasible. but attempting this approach is worthwhile. Monitoring blood 
volume changes with water absorption is not critical to the success of our imaging 
approach, but it has the potential to make the overall imaging approach more powerful if 
it works. 

We will increase the power available from our LED array. The images in the 
Preliminary Studies section were obtained using 21 LEDs at each wavelength. We will 
increase the power by a factor of 20 by using more LEDs and increasing their brightness 
by operating them at higher drive currents. We have shown in burn-in tests that the 
LEDs can be operated significantly above their typical operating currents for many 
weeks without incurring problems. The LEDs are turned on for only short periods during 
imaging at each wavelength, thereby increasing the practicality of higher current 
operation without LED damage. 

We will substitute a more sensitive camera for the one cited in the Preliminary 
Studies section. We will purchase a Retiga 1350 EX camera produced by Q-lmaging. 
This CCD camera is approximately 2 times more sensitive in the NIR than our current 
camera. In addition, the camera-sensitive area is 4 times larger. These two 
improvements will lead to an overall enhancement in camera sensitivity of roughly a 
factor of 8. In combination, the increased illumination and more sensitive camera should 
improve overall system sensitivity by more than 100 times. 

Aim 2. Perform Human Studies 

We will perform optical imaging on human subjects to test the approach. These 
studies will be perfomned at the Stanford University Medical Center (approximately 2 
miles from SRI) under the overall supervision of Dr. Birdwell. Associate Professor of 
Radiology. Dr. Birdweirs assistant will be responsible for recruiting volunteers, 
maintaining records, and coordinating scheduling. The imaging system will be operated 
by a postdoctoral associate from SRI. To provide the patients being imaged with a 
greater sense of ease with the imaging process, the postdoctoral associate will be a 
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woman. This imaging may be performed in either the stereotactic biopsy room or the 
ultrasound-guided biopsy room at Stanford. 

Subjects will be recruited from women referred for biopsy, either on the basis of 
screening or on the presence of a suspicious lump. Optical imaging will be performed 
before the biopsy to avoid any Influence the biopsy procedure might have on imaging 
measurement and interpretation. The imaging will be performed using only one or two 
inhalation protocols so that the total imaging takes only a few minutes. Imaging will 
cease if the patient feels uncomfortable or wishes to stop for any other reason. Imaging 
will be performed on 10 subjects in the first year and 15 patients in the second and third 
years, for a total of 40 patients. If the sensitivity and specificity are high, measurements 
on this number will enable us to calculate 95% confidence Intervals for specificity and 
sensitivity with a half-width of approximately 15%. Volunteers for the imaging will be 
given a remuneration of $50. 

Different imaging methods may be used for differential vasoactive imaging of the 
breast. The imaging may be performed with or without compression and with or without 
immersion. We hypothesize that optimal imaging will entail using at least mild 
compression and immersion. Compression would be advantageous for two reasons: 
first, with compression the total imaging distance is less, leading to a higher S/N. and 
hence increasing the likelihood of detecting a smaller tumor. Second, X-ray 
mammography uses compression. We believe that combining optical imaging with X-ray 
imaging would provide a powerful method— given their low-cost and the possibility that 
both imaging techniques could be perfomned simultaneously. In fact, with the transition 
to digital mammography using semiconductor-based cameras, the imaging could share 
the same detector. That combination would lead to an improvement in sensitivity and 
specificity over either modality alone. Achieving this improvement would require 
coregtstration of images from the two modalities, which could be achieved most 
practically if compression were used for each imaging method. 

From our experience with the animal imaging described in the Preliminary Studies 
section and other in vivo imaging,^ we believe that the sensitivity of the imaging will be 
highest when immersion is used. With immersion, all portions of the breast are imaged, 
with nearly the same illumination reaching the detector and providing more optimal use 
of the dynamic range of the camera. That is, the entire image may be acquired with a 
high level of illumination, and hence high S/N. For the nonimmersed breast, variations in 
the transmitted light intensity across the breast will be large. To avoid camera saturation 
in the thinnest regions, low light levels will be obtained in the thicker regions. Thus, the 
thicker regions will have a lower S/N, and worse imaging results. Researchers have 
used the phase measurement available with frequency domain measurements to 
perform correction for edge effects.®^ Immersion achieves a similar goal. 

Immersion can be achieved in two ways: first, the patient lies prone on a table similar 
to a stereotactic breast biopsy table with the breast immersed in a matching medium 
below.®-®^ A second approach is to surround the breast with a doughnut-shaped 
transparent bag containing a tissue phantom liquid. The bag would be filled to a slight 
overpressure to press against the breast in a manner similar to a blood pressure cuff, 
except that the overpressure would be much less. This method would achieve the same 
advantage of immersion but with less preparation and cleanup required. We plan to use 
the second Immersion method, if we find significant problems, we will use the method 
with the breast pendant in the matching medium. For immersion-based imaging, a new 
bag with fresh immersion medium will be used for each patient. The immersion medium 
will be maintained at 37''C. 
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Imaging will be performed only on patients who already have suspicious lesions. 
However, optical imaging could identify another suspicious breast region that is not 
identified by the imaging modality being used for the biopsy (X-ray or ultrasound). 
Because the diagnostic capability of the optical imaging is still unproven, further imaging 
of a lesion identified only by optical imaging will not be covered by the patient's medical 
insurance. To address this eventuality, the budget includes the cost for MRI for follow- 
up on lesions seen only with optical imaging. If MRI does not detect the lesion that has 
been detected optically, the patient will be scheduled for a follow-up examination after 6 
months. 

The relative sensitivity and specificity of a diagnostic method depend on the criteria 
used. Relevant criteria Include percentage change in hemoglobin content and 
hemoglobin oxygenation, and the relative signs (i.e., did each increase or decrease). As 
with any diagnostic method, by varying the criteria used either sensitivity or specificity 
can be made high, but at the expense of the other dimension. To assist in the analysis 
of the data, we will use a receiver operating characteristic (ROC) curve, which plots 
sensitivity versus false positive fraction; the free parameter is the criterion or threshold 
used for diagnosis. The area under the ROC curve gives a measure of the quality of the 
method; an area near 1 is desirable. We will prepare ROC curves for each contrast 
mechanism and for the contrast mechanisms in conjunction. 

To test the null hypothesis of no association between optical mammography results 
and biopsy status, we will perform a sign test, wherein a correct prediction (i.e.. a patient 
with agreement between the optical mammography and biopsy results) is coded as a 1 
and an incorrect prediction is coded as a 0. Under the null hypothesis (and conditional 
on observing 8 patients with cancers as determined by biopsies), the number of correct 
predictions follows a distribution that is the combination of a Binomial(8. 0.2) and a 
Binomial(32. 0.8) distribution (i.e., 8 flips of a coin with a probability of 20% of being 
correct and 32 flips of a coin with a probability of 80% of being correct). Under the 
conservative alternative hypothesis, the number of correct predictions follows a 
distribution that is the combination of a Binomial(8, 0.8) + Binomial(32. 0.8) = 
Binomial(40, 0.8) distribution (that is. 40 flips of a coin with a probability of 80% of being 
correct). 

We reject the null hypothesis if we obtain more correct predictions than we would 
expect under the null hypothesis (i.e., the number of con^ect predictions exceeds the 95- 
th percentile of the distribution under the null hypothesis). Monte Carlo simulation (1000 
trials) shows that the 95-th percentile of the distribution under the null hypothesis is 32 
(i.e., there is only a 4.3% probability that under the null hypothesis there will be 32 or 
more connect predictions). Consequently we reject the null hypothesis if we find 32 or 
more correct predictions. Monte Carlo simulation (1000 trials) under the conservative 
alternative hypothesis shows that there is a 60% probability of observing 32 or more 
correct predictions. Therefore, the power is 60% under the conservative hypothesis. If 
the results in humans are as good as the results in animals, then our specificity and 
sensitivity will each be approximately 90%. Under this alternative hypothesis, the power 
of the study is approximately 98% (i.e.. there is a 98% probability that we will find 32 or 
more correct predictions out of 40 patients). 

Aim 3. Perform Animal Model Studies 

We will supplement our human studies with limited animal model studies to achieve 
better understanding of findings for humans, of differences between animals and 
humans, and of variation between tumor types. Work with animals will also help 
minimize inconvenience to subjects and will allow greater flexibility. For example, we 
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can perform studies over a longer period, perform more measurements in a single day, 
more readily repeat measurements on different days, compare results with identical 
cancer types in different animals, compare results on different cancer types, and 
measure small tumor sizes. 

We will perform imaging in a number of additional animal cancer models to improve 
our understanding of the range of response possible with vasoactive differential imaging. 
Understanding how the response differs among tumor types will enable us to adapt the 
image interpretation appropriately. Ideally, each different tumor type or stage of tumor 
growth will display sufficiently reliable characteristics to allow tentative identification of 
tumor type or assist In cancer staging. The animal model will help us reproduce results 
found in the human studies and to vary those results. 

Animal imaging will also allow us to refine our instrumentation and data analysis. We 
can perform measurements on the same animal at different times during tumor 
development or use different gas inhalation protocols. This wori< will enhance our ability 
to understand how the imaging signatures change during tumor development and how 
different gas inhalation protocols change the imaging response for the same tumor and 
animal. 

We have experience with a number of breast cancer models for this animal model 
work. To gain a greater sense of the variability of the differential vasoactive response, 
we will compare an orthotopic model with various ectopic tumors (both syngeneic and 
xenograft) and orthotopic tumors In lactating rats. Comparison of the rat orthotopic 
model with measurements on ectopic mouse tumors will provide a good sense of the 
variability of the differential vasoactive response. We will perform contrast 
measurements on three rat models (orthotopic mammary tumors, orthotopic mammary 
tumors in lactating mammary glands, and ectopic mammary tumors) and three mouse 
models (4T1, MCF*7. and MDA-231). We will use a total of 30 rats and 30 mice to 
obtain modest confidence Intervals for our preliminary studies (about 15%). 

Rat Orthotopic Breast Cancer Model 

To construct a cancer tumor model that gives representative results with improved 
reliabiliW, we will use orthotopic tumors with the syngeneic mammary adenocarcinoma 
R3230. '^^ Orthotopic tumors are considered to be more similar to human mammary 
carcinomas In regard to histology and structure than ectopic tumors t>ecause they are 
grown in the site of origin (i.e., the mammary gland). In addition, the makeup of the 
mammary gland — a collection of inactive lobular glands with abundant stroma and fat 
tissue — is quite similar across different species. Accordingly, we will use syngeneic 
orthotopic mammary adenocarcinoma R3230 in female Fisher 344 rats. We will inject 
5 million cells into the subcutaneous tissue of the thigh. Once the donor tumor reaches 
a volume of 400--600 mm^, the tumor will be surgically removed. We will surgically 
implant 1*mm^ tumor pieces from the donor tumor into the mammary fat pad of recipient 
female F344 rats. We will use one distal mammary gland on the mammary line and use 
the other glands as controls. 

Orthotopic Mammary Tumors in Lactating Mammary Gland 
We will use lactating mammary glands as one of the sources of normal tissue for 
comparison of the tumor results. The multifold enhancement of the blood supply in the 
lactating mammary gland will challenge the ability of the imaging approach to distinguish 
normal and cancerous tissue. Lactation signifies the most active state of the mammary 
gland because the blood supply is enhanced multiple times to supply functional glands 
that are secreting milk. As a result, imaging is much more difficult because of the 
extensive noise induced by that activity inside the mammary gland, especially if a 
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growing tumor is present. We will use pregnant female Fisher rats for these 
experiments. With surgical implantation of 1-mm^ tumor pieces from the donor tumor, 
tumor growth will proceed at a rapid rate. This is important because lactation will 
continue for no more than 3 to 4 weeks postpartum. Only one gland will be selected for 
the tumor injection, with the other glands used as controls. 

Ectopic Tumors 

We will also study ectopic tumors using the same cancer cell line. Study findings will 
provide information on the relative roles of the cancer cell line and of the tumor site on 
the nature of the vasoactive imaging. For the ectopic model, we will surgically implant 
tumor pieces approximately 1 mm^ in size subcutaneously into the thighs of Fisher 344 
rats. 

Mouse Mammary Carcinoma 4T1 Ceii Line 

The 4T1 cell line^^^'^^® is a very aggressive metastatic mouse mammary cancer cell 
line. The tumors will be generated through cell injection of 1 million cells in 
subcutaneous tissue of the thigh region of the female Balb/c mouse. The tumor volume 
will be measured twice weekly. The imaging will be carried out in tumors having 
volumes ranging from 200-800 mm^. 

Human Cancer Cell Unes-Xenografts 

We will use two commonly used human breast cancer cell lines, MCF-7^^®-^^^ and MDA- 
231.^^® to generate orthotopic (mammary fat pad) tumors in nude mice. Since MCF-7 is an 
estrogen-dependent cell line, the animals will also be injected with an estrogen pellet over 
the dorsum. The tumor volume will be measured twice weekly. The imaging will be carried 
out in tumors having volumes ranging from 200-800 mm^. 

Animal Care and Handling 

The SRI animal research facility is AAALAC-accredited. The temperature is 
maintained at 72 ± 5^F, with relative humidity at 35 — ^70% and a 12-h light/dark cycle. 
Rodents have certified Purina rodent chow and drinking water freely available. The 
source water is recirculated, deionized, UV-treated, and 5-^m filtered. It is analyzed 
quarterly for metals and trace contaminants. 

Each rat will be sacrificed following imaging and before undue pain from the tumor 
occurs. Rodents will be sacrificed by AVMA-approved euthanasia solutions (i.e.. sodium 
pentobarbital [200 mg/kg LP.] or Buthanasia [150 mg/kg LP.]). If tumor volumes reach > 
1 ,000 mm^ or if weight loss exceeds 20%, rats will be sacrificed regardless of whether 
the planned imaging is complete. 

Animal Imaging 

Animals will be anesthetized with sodium pentobarbital (40 mg/kg LP.) before 
imaging to minimize trauma and to facilitate performance of the imaging. Imaging will be 
perfomied using a method similar to that used for the preliminary studies. The 
anesthetized animals will be taped to a plastic paddle and immersed in the tissue 
phantom, which will be maintained at 37^0. It is known that anesthetization can cause 
physiological changes^^ that might affect our Imaging measurements; this is especially 
the case for the ketamine/xylazine combination, which we have not opted to use. 
However, changes associated with sodium pentobarbital are expected to be fewer than 
the physiological changes associated with an animal under stress from constraint. 

The animals will be monitored regularly and the tumor volume measured twice 
weekly. We will carry out optical imaging at three time points corresponding to different 
tumor sizes (200 mm^, 500 mm^ and 1000 mm^). The animals will be anesthetized with 
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35-45 mg/kg of pentobarbital before the start of the imaging procedure. After the imaging 
session, the tumors will be removed and fixed in 10% fonmalin. Histological analyses will 
be carried out on H&E (hematoxylin and eosin) stained slides. We will conduct 
semiquantitative observations of certain well-defined morphological changes that occur 
as the mammary carcinoma grows, such as necrosis, cystic changes, and development 
of solid areas, and compare them with the optical imaging results. Such changes should 
be observable optically. Necrosis and cystic changes should both have low hemoglobin 
levels, compared with an actively growing tumor, and so be clearly identifiable. 

* Data Analysis 

The measured properties for tumors will be compared with those of neighboring 
tissue and lactating mammary glands in the animals. To verify the anatomy observed in 
the images, tissues from rats will be embedded in paraffin and processed for 
pathological examination. We will carry out hematoxylin and eosin stain for general 
histology and immunohistochemistry for blood vessel identification to gain an 
understanding of the results from optical imaging and actual tumor morphology at the 
last imaging session. The nature of the tumor vasculature from the pathology will be 
compared with the imaging results to determine how it contributes to the images 
obtained. This information will be helpful in continuing development of this imaging 
approach. Drs. Amin, Haroon, and Dewhirst will provide advice on cancer biology 
relevant to the images. Dr. Amin is at SRI and will be performing the animal model work. 
Dr. Haroon, a consultant for this project, has experience with the proposed animal 
models and the role of hypoxia in cancer. His offices are near SRI, and he is readily 
available for consultations. Dr. Dewhirst, also a consultant, will provide additional 
expertise on the role of O2 transport in tumors. He is performing related research using 
MRI instead of optical imaging. 

Aim 4. Study Gas Protocols 

We will examine other protocols for gas inhalation. Measurements will be performed 
with varying inhalation gas composition and administration time. Gas mixtures will be 
produced on demand using computer-controlled gas flow controllers. Three gases will 
be used to produce these mixtures: nitrogen, O2. and CO2. In this way we can rapidly 
alternate among gas compositions, continuously varying the levels of CO2 and O2 in a 
nitrogen buffer, or create carbogen. These studies are of interest because CO2 and O2 
have opposing effects on vasculature (vasodilation versus vasoconstriction, 
respectively). Using these two mechanisms in opposition or in alternation should 
produce potentially useful results from the differential vasoactive imaging. For example, 
elevated CO2 levels may be administered for a period of a minute, followed rapidly by a 
period of elevated O2. The same protocol could be repeated with a small overlap 
between the elevated CO2 and O2 levels. With the computer-controlled flow controllers, 
we can sequentially administer different gas mixtures to the same individual, taking care 
that the vasculature recovers sufficiently between the changes. We expect that these 
studies will result in more effective discrimination between cancerous and noncancerous 
tissue, and perhaps provide means of distinguishing among different tumor types. 

We will perform studies of gas inhalation protocols on both animals and humans. 
With animals, we can perform studies over a longer period, conduct more 
measurements in a single day, more readily repeat measurements on different days, 
compare results with identical cancer types in different animals, and measure small 
tumor sizes. The protocols are likely to be different for humans because of the different 
respiratory rate, heart rate, size, and the fact that the animals are anesthetized and the 
humans are not. 
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Animal studies will be performed on different animal cancer cell lines to determine 
whether we can distinguish among those lines. If we find certain protocols accentuate 
the difference among cell lines or animals, we will modify that protocol to see whether 
the degree of discrimination can be increased. If certain protocols reduce discrimination, 
we will modify those protocols to determine whether the degree of discrimination can be 
decreased while still maintaining the discrimination between cancerous and 
noncancerous tissue. Protocols with and without strong discrimination are both likely to 
be of clinical value. Protocols that do not discriminate strongly among tumor types may 
prove more useful for imaging related to screening, when the primary interest is cancer 
detection. Protocols with stronger discrimination among different tumor types might be 
used when suspicious lesions are detected. The more promising gas protocols found for 
the animals will be used in the clinical studies, provided that the gas compositions 
remain within safe levels for humans. 

Aim 5. Develop Image Analysis Tools 

The measurements acquired for differential vasoactive imaging comprise three- 
dimensional (3-D) datasets as illustrated in Figure 1 1 . The two spatial dimensions and 
one temporal dimension differ from other 3-D imaging modalities such as MRI or 
computed tomography (CT), which have three spatial dimensions. For those imaging 
modalities, visualization tools often create 2-D images as cross sections through the 3-D 
data set. Regions of interest can be probed by changing the orientation of the cross 
section. This same approach is also used for ultrasound imaging, where the ultrasound 
technician achieves different sections in real time by changing the orientation of the 
ultrasound probe. 

Visualization of data such as in Figure 1 1 can be performed by taking cross sections 
at different orientations. In fact Figure 2 and Figures 3 and 4 are examples of a cross 
section and a line section through such a data set at constant time and position, 
respectively. However, both the spatial pattern (such as Figure 2b) and the temporal 
pattem (such as Figure 3) are necessary to define features in this data set. 
Simultaneously capturing both of these features requires a different sort of image 
analysis tool. One such tool is PCA. 




Figure 1 1 . Form of 3-D data for differential vasoactive imaging. 

We will apply PCA similar to that described in the Preliminary Studies section. We 
will adapt the approach to allow automated data processing of temporal image data sets 
of oxyhemoglobin, deoxyhemoglobin, and total hemoglobin, and change in O2 content. 
To improve the results, we will study methods such as spatial and temporal averaging 
conditioned on the image features and the use of a priori information such as the 
temporal profile of the gas inhalation protocol. For breast imaging, we expect the 
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number of eigen images may be larger. We will study methods for classifying the eigen 
images (e.g., by tumor type, other feature such as blood vessels). 

Future Work 

Following this TRC project we will have a very good indication of whether DVOI is 
effective for breast cancer detection. If the technique shows good sensitivity and 
specificity, we plan to continue clinical development of the approach. Depending on the 
outcome of the human studies, this future work may involve further imaging system 
refinements, human studies with larger numbers of subjects, human studies targeted at 
specific imaging needs such as screening young women with known propensity for 
developing breast cancer, or development of a combined imaging modality such as 
optical and x-ray imaging. 
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